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A fractionation model based on three lognormal particle size distributions 
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In this paper, a new model is proposed to calculate distribution of fission products in particles of different 
sizes. The model sensitivity to the effective volume and mass of vaporized soil particles is examined. Compared 
with other fractionation models, the new method has a much better performance in calculating rg9,95, but the 
calculated cumulative activity fraction for particles in diameters over 100 um is in between the results using the 
F-T and G-X models. It is concluded that in a near surface nuclear explosion radioactivity is mainly distributed 
in soil particles which have not been vaporized, and according to the Henry’s law and ideal gas law, rg9,95 may 
vary in larger particles when effective volume of the fireball is changed. 
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I. INTRODUCTION 


A nuclear explosion generally produces five effects: blast 
and shock, thermal radiation, electromagnetic radiation, ini- 
tial nuclear radiation, and residual nuclear radiation (fallout) 
[1]. Fallout plays an important role especially in ground nu- 
clear explosion by creating a concentrated radioactive hazard 
causing potential damages to ecosystem in a relatively long 
period of time. Radioactive fallout must be evaluated to limit 
the ionization damages and provide guidance for formulating 
emergency plans in response to nuclear explosions. The to- 
tal amount of radioisotopes and their distribution in particles 
of different sizes are determined by parameters of nuclear 
explosion, such as weapon yield, explosion height and soil 
type of ground zero. Therefore, computer simulation of the 
formation of radioactive particles in various explosions and 
quantitative analysis of the radioisotope distribution in differ- 
ent particles play an important role in fallout forecasting. It 
was found that the influence coefficient of fractionation in ra- 
dioactive fallout is over 2 according to researches conducted 
in Air Force Institute of Technology (AFIT), but fractionation 
models are still important for better understanding of explo- 
sion phenomena and evaluation of validity of samples [2]. 

Freiling-Tompkins fractionation model (also known as 
Modified Radial Power Law), Miller’s thermodynamic equi- 
librium model and Martin’s diffusion based G-X model are 
typical fractionation models [3-5]. The G-X model was de- 
signed by Martin in 1980s, in which fallout particles are di- 
vided into two distinct types: glassy and crystalline. How- 
ever, the total specific activity of glassy and crystalline parti- 
cles calculated with the G-X model varies much more greatly 
than the results of test observation. This can be attributed 
mainly to the different particle position in the fireball. All 
the particles are drawn into fireball at an early time, and dif- 
ferent locations in the fireball decide the specific activities 
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because of rapid changes in temperature and fission product 
concentration in the fireball. Also, some fine particles may 
come into being from condensation of the vaporized soil on 
separate atoms/molecules. After the soil particles begin to so- 
lidify, soil particles of smaller diameters are still drawn into 
cloud by after-winds, and are entirely un-melted. 

Considering the particle formation process, the fallout par- 
ticles can be divided into three types: 1) particles from the de- 
vice debris and vaporized soil, 2) particles representing com- 
pletely or partially melted soil, and 3) particles representing 
soil never melted. Based on the particle size and formation 
process, the three types of particles are referred to as fine par- 
ticles, glassy particles and crystalline particles, respectively. 
With this classification, a new fractionation model was con- 
structed to describe the distribution of radioactive nuclides 
in these particles. Taking into account that effective volume 
and the mass fraction of fine particles cannot be determined 
clearly, influence of the result from effective volume and mass 
of fine particles was also examined. 


Il. METHOD AND PARAMETER SELECTION 


A. Method 


Calculation proceeds on one decay chain at a time when 
the fireball temperature is about 3500 K. After this point the 
isotopes are assumed to condense on surface of the particles 
according to Henry’s law and diffuse into the fallout particles. 
The isotopes decay with the time and the temperature falls. 
The process is repeated until the temperature is close to the 
solidification temperature for the glass. After this point, the 
fission products in gaseous state will be distributed on sur- 
face of all types of particles. The amounts of the nuclides 
distributed in particles of different sizes and types are propor- 
tional to the total surface area of the particle group. 

In the first stage, only agglomeration between fine particles 
should be considered. Since all the fine particles are small 
enough to be considered nuclides volume-distributed, if we 
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ignore the diffusion process in the fine particles, there is no 
relation between the distributions of fission products and the 
fine particle growth process. To simplify the model, we sup- 
pose that the fine particles of a constant mass at the first stage 
grow up mainly by agglomeration. The diffusion equation for 
spheres with only radial diffusion is 


Oc; 07 ci 2 Oc; 
D; (Zs + 2 ; (1) 


ot 


where c; is the concentration of fission product 7, t is time of 
the diffusion process, r is the radial coordinate, and D; is the 
diffusion constant. According to the Henry’s law equation, 
the partial pressure P; of fission product 2 is 


where k; is the Henry Constant and x; is concentration of 
fission product 2 at the particle surface. 

By dividing the first stage into many small steps, at Step 1 
we have the following mass balance by combing the diffusion 
equation, Henry’s law and ideal gas law: 
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where subscript ¿ denotes the fission product; Y; is its mass 
i in the first time period; c; is its surface layer concentration 
in the soil and k; is its Henry’s Law constant, at temperature 
of Step 1; R is its gas constant; M; is its mole weight; Mms 
is the mass of soil in particle size class s; p the number of 
particle size classes; m, the mass of fine particles; and tı is 
the first time period. On right side of Eq. (3), the first term 
is the amount of fission product 7 in gas phase, the second 
term is the amount of fission product 7 diffused into melted 
soil particles by solving diffusion equation, and the third term 
is the amount of fission product 2 diffused into fine particles 
formed mainly by vaporized soil. 

The immediate precursor of isotope ¿ which diffused into 
particles in time interval At; may have decayed into isotope 
i. The influence to diffusion of isotope 2 in Step 2 by the men- 
tioned process can be ignored if the time step is small enough. 
The mass balance equation after the second time period is 
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where Y; is the mass of fission product 7 and c;2 is the surface 
layer concentration of i in the soil, in Step 2; A; is the decay 
constant of i; and A;_, is the decay constant of precursor of 
isotope 7. 

This calculation goes iteratively until the temperature is 
close to the solidification temperature for the glass. In the 
second stage, the distribution of fission product 2 in particles 
with different sizes and types is given by 


Ymsi = Ymy $ Sur fj i Ía (7) 
myv; 3 exp { (—2.5- In?(a;)) } 
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j= 
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where Yms; is the mass of isotope 7 on the surface of the sth 
size and j" type particles, Y,,, the mass of gaseous istope i 
after stage 1, Sur f; the total surface fraction of the j® type, 
m; the total mass of the j t type particles, vj the mean volume 
of the j"" type particles and gj the standard deviation. 


B. Particle size distribution 


Heft suggested that nuclear debris might be described by 
a linear combination of several lognormal distributions [6], 
which are distributions for each type of particles generated 
by the burst. The subsurface mass distribution is given by 


10 = 32 (Sie) om e] 


where g; is the percentage of mass for the i" type of particles, 
o; is the standard deviation, and ao; is the natural logarithm 
of average diameter. 

Based on Delfic’s lognormal distribution and Baker’s air- 
borne particle distribution, and with the knowledge that soil 
particles near ground zero of the Small Boy atom bomb are 
below 150 um in diameter, the values of o; and ao; in this 
research are given in Table 1. 


TABLE 1. The values of g; and ao; in this research 


Particle types Qoi Gi 
Melted soil particles 0.407 4.0 
Un-melted soil particles 4.000 2.0 
Vaporized soil particles 0.200 2.0 


C. Parameters selection 


The main fireball parameters in this research are the vol- 
ume and temperature during the fireball cooling to the solid- 
ification temperature of glass. The temperature is selected 
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from the Small Boy data. The fireball size varied little af- 
ter the fireball breakaway, so the averaged fireball volume is 
10° m?/kt in this research. 

Norman et al. did many studies to determine Henry’s Law 
constants for those atomic species of interest in fallout re- 
search [7]. Norman and others contributed a lot in determin- 
ing diffusion constants of various soil types. The soil type of 
interest here is the CaOQ—Al,03,—SiO, soil, which is repre- 
sentative of common soils. The values of Henry’s Law con- 
stants and diffusion coefficients used in this research are listed 
in Ref. [3]. 

Fission yield data and radioactive decay data are required 
for fission product dynamic calculation. Independent fission 
yield data together with radioactive decay data were selected 
from ENDF BVII.1 [8]. 


Ill. RESULTS 
A. Fractionation results 


In a ground explosion, about 5000 t of earth per kt of yield 
are released into the atmosphere and, of that quantity, 180 ~ 
200 t are fused. The vaporized soil is 1.5 ~ 25t per kt. As the 
Small Boy is a nuclear explosion of low yield, 3t vaporized 
soil particles and 1 kt crystalline particles are assumed and 
200 t melted soil particles is determined [9]. 

To evaluate the model, the plots of rgg95 as a function of 
particle size is useful. The ratio r; j is defined as 


rig = fil fi, (10) 


fi =ai/Yi, (1) 


where a; is the number of atoms of nuclide 7 in the sample 
corresponding to one fission and Y; is the yield of nuclide 7 
per fission. Log-log plots of r; 95 versus rgg95 values can be 
prepared and fitted with linear least squares to determine the 
correlation slopes. 

The rgg.95 and the correlation slopes calculated by this 
model and G-X model are compared to the Small Boy data 
in Fig. | and Table 2. The cumulate activity-size distribution 
calculated by this model with 3t vaporized soil particles and 
1 kt crystalline particles is compared in Fig. 2 with the results 
using G-X and F-T models. 

Figure 1 shows that our method better matches the Small 
Boy data above 100um. The Small Boy data are mainly 
above 20 um, so the new method performs much better in cal- 
culating rgg.95. Fig. 2 shows that our results of the cumulate 
activity fraction of particles over 100 um diameters, which 
mainly grounded in an area near ground zero after explosion, 
are between the results of the F-T and G-X models. But the 
logarithmic correlation slopes calculated with both models 
shows poor consistency with the Small Boy data, but the new 
model shows a poorer consistency in logarithmic correlation 
slopes of 131. Log-log plots of 7131.95 versus rgg9 95 of G-X 
model in Ref. [3] are cited, and the linear fitting gives a slope 
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Fig. 1. (Color online) rgg,95 vs. particle size for different methods. 
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Fig. 2. (Color online) Cumulate activity-size distribution for differ- 
ent models. 


TABLE 2. The values of g; and ao; in this research 


Atomic mass Small Boy G-X This work 
89 1.00 1.00 1.00 
90 0.73 0.83 0.73 
91 0.51 0.64 0.39 
95 0.00 0.00 0.00 
99 0.04 -0.07 0.09 
131 0.84 0.58 0.36 
132 0.90 0.61 0.67 
136 0.65 1.39 1.35 
137 1.19 0.91 1.02 
140 0.58 0.34 0.34 
144 0.03 -0.14 0.01 


of about 0.40 which varies much with the value given in Ta- 
ble 2 and is nearly equal with result of this model. As we can- 
not find more information about how Martin fitted the plots 
and logarithmic correlation slopes vary little when we change 
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effective volume and fine particle mass fraction, we do not 
pay more attention to the logarithmic correlation slopes in 
later research. The current deviation of slopes of mass chain 
91 and 140, where elements mainly affecting fractionation in 
the chain are mixed behavior elements whose distribution in 
fallout particles is appreciable gradient or some bulk loading, 
may suggest that diffusion constants used are overestimated. 


B. Sensitivity analysis 
1. Model sensitivity to effective volume 


Since the fireball experiences a rapid temperature and con- 
centration fall, the effective volume in which the condensa- 
tion process is involved is smaller than the physical volume 
of the fireball. The rgg 95 and activity-size distribution were 
shown in Figs. 3 and 4, with 20% and 100% as the fraction of 
effective volume in the fireball. 


10? 


Fraction of effective volume 
—e— 20% 


—=— 100% 


10° 


10-5 1074 10 


Diameter (m) 


Fig. 3. (Color online) Dependency of particle size on rgg,95. 
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Fig. 4. (Color online) Cumulate activity-size distribution. 


For particles of diameters over 1 mm, as shown in Fig. 3, 
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Fig. 5. (Color online) rs9,95 vs particle size for 0-15 t mass of va- 
porized soil particles. 


1.0 4 
Vaporized soil particle mass ra 

0.8 xy 
=i- 0t =t = Z 


=: P 
A a 


Cumulate activity fraction 
© 
N 


2 
ro) 


103 10+ 10° 
Low limit of diameter (m) 


Fig. 6. (Color online) Cumulate activity-size distribution for 0-15 t 
mass of vaporized soil particles. 


the rgg,95 differ from each other for 20% and 100% fraction 
of effective volume, while the cumulate activity-size distribu- 
tions remain unchanged. According to Henry’s law and ideal 
gas law, the fraction of fission products condensed in parti- 
cles or still in gaseous state changes when effective volume 
varies. 


2. Model sensitivity to mass of vaporized soil particles 


In surface nuclear explosions, fine particles may form in 
the process of condensation on separate atoms/molecules. 
But mass fraction of the particles cannot be determined due to 
uncertainty of vaporized soil mass and the condensation pro- 
cess competition between the fine particles and melted soil 
particles. About 40t fine particles may form in an explosion 
of 1.65 kt weapon yield, but just 15 t fine particles was used as 
a maximum value because of the low yield and about 3 m ex- 
plosion height. The rgg,95 and activity-size distribution were 
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calculated for vaporized soil particle masses of Ot, 3t and 
15t (Figs. 5 and 6). The results show that the mass of vapor- 
ized soil particles has no effect on the rgg 95—size relation and 
the cumulate activity—size distribution hardly changes with 
the vaporized soil particle mass, as melted and un-melted soil 
particles involved are much more than the vaporized soil par- 
ticles. 


IV. CONCLUSION 


Large quantity of soil entrained into the fireball by blast 
and updraft winds is not melted. We can assume that the un- 
melted soil particles do not play an important role in the con- 
densation process of fission product on the particles, but they 
offer large surface areas for the gaseous fission product after 
condensation of the melted soil. Based on the Small Boy data 
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